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The MgxZni _xO alloy system may provide an optically tunable family of wide band gap materials that can be 
used in various UV luminescences, absorption, lighting, and display applications. A systematic investigation 
of the MgO-ZnO system using ab initio evolutionary simulations shows that Mg^Zni-^O alloys exist in 
ordered ground-state structures at pressures above about 6.5 GPa. Detailed enthalpy calculations for the 
most stable structures allowed us to construct the pressure-composition phase diagram. In the entire 
composition, no phase transition from wurzite to rock-salt takes place with increasing Mg content. We also 
found two different slops occur at near x = 0.75 of curves for different pressures, and the band gaps of 
high pressure ground-state MgxZni_xO alloys at the Mg concentration of x > 0.75 increase more rapidly 
than X < 0.75. 

To realize high-performance ZnO -based optoelectronic devices, two important requirements are necessary: 
one isp-type doping of ZnO, and the other is modulation of the band gap (£g)^'^. ZnO can alloy with MgO to 
form the ternary compound MgxZni_xO to extend the energy band gap, and therefore the detection 
spectrum into the shorter wavelength region. While p-type doping of ZnO is studied intensively, the latter has 
been demonstrated by the development research on MgxZni_xO allowing modulation of band gap in a wide 
range, from 3.34 to 7.8 eV^. Although the topic of ZnO has been extensively researched, less is known concerning 
the detailed structures of the MgxZni_xO alloy system. 

Fabrication and characterization of MgxZni_xO alloy are important from the viewpoint of band gap modu- 
lation as well as oip-n junction. MgxZni_xO thin films and nanostructures that include nanocrystals as well as 
core-shell structures have recently been studied with the objective of achieving a viable alloy family with tunable 
band gap and luminescence at the UV range'*"^\ Sans et al.^ investigated the dependence of the phase transition on 
the hydrostatic pressure and found that the wurtzite to rock-salt transition pressure is observed to decrease from 
9.5 ± 0.2 pure (ZnO) to 7.0 ± 0.2 GPa (for x = 0.13), with an almost linear dependence on the Mg content. 
Difficulties for the growth of high quality MgxZni_xO films are partly from the fact that MgO (cubic) and ZnO 
(hexagonal) have different crystal structures at normal conditions. 

In the last decade, there are numerous theoretical studies on the MgO -ZnO alloys^^"^°. Sanati et aU^ showed that 
isostructural MgO-ZnO alloys are stable under certain conditions using the cluster expansion (CE) method, but 
they also concluded that if MgO and ZnO can adopt their own crystal structures {Bl and M, respectively), the alloy 
is predicted to phase separate. This means that MgxZni_xO alloys are not thermodynamically stable, consistent 
with a rather low observed solid solubility limit for Mg in ZnO. The cluster-expansion method was also used by 
Seko et aV^ for investigating phase transitions, including vibrational effects through lattice dynamics calculations. 
Very recently, combining the CASP-CE with a systematic set of first -principles total energies, Liu et a/.^^ exhaus- 
tively searched for the ground-state structures of MgO-ZnO alloys, and found a few structures as yet unreported. 
Several groups^^"^^ reported their theoretical studies in which the random MgxZni_xO alloy have been simulated 
using special quasirandom structures (SQS) approach^^'^^. There exist other theoretical studies^°"^^ in which 
alloying of ZnO and MgO proceeds by substituting Mg atoms by Zn atoms in the cubic rock-salt structure or 
vice versa in the hexagonal wurtzite structure, and also exist several theoretical reports^^"^^ on MgZnOa in various 
model systems. An extensive systematic study of structural properties of MgxZni_xO was performed by Fan et aV^ 
using a supercell approach within the local-density approximation. In another theoretical study^°, alloying of MgO 
and ZnO is described within coherent-potential approximation (CPA)'*^ 

Bulk and undoped ZnO prefers the hexagonal wurtzite (54) structure under normal conditions and transforms 
into a cubic rocksalt {Bl) structure at a pressure in the vicinity of 9 GPa reported by some experiments^^"^^, while 
MgO adopts the cubic rocksalt structure up to the highest pressure of 227 GPa^^. Generally, a large crystal 
structure difference between the wurtzite-hexagonal ZnO (M) and the rock-salt-cubic MgO {Bl) can cause 



SCIENTIFIC REPORT; | 4 : 5759 | DOI: 1 0. 1 038/srep05759 



1 



B3-ZnO+Bl-MgO 
B4-ZnO+Bl-MgO 
^Bl-ZnO+Bl-MgO 




0.2 0.4 0.6 0.8 

X |MgO/(MgO+ZnO) ratio) 



B4-ZnO+Bl-IVIgO 



-MgZn304 
-MgZn02 
- Mg3Zn04 




7.3 GPa 



(ci) 
0.20 
0.15 
0.10 
0.05 
0.00 

-0.05 



0.4 0.6 0.8 

X [MgO/(MgO+ZnO) ratio] 









-■-MgZn304 


B3-ZnO+Bl-MgO 


-#- MgZn02 




-A- Mg3Zn04 




7.5 GPa 






5.9 GPa 





Pressure (GPa) 



Pressure (GPa) 



Figure 1 | Thermal stability of MgO-ZnO alloys. a,b Predicted formation enthalpy of MgO-ZnO alloys with respect to decomposition into Bl-MgO and 
ZnO (B4, B3, or Bl) at (a) 0 GPa and (b) 8 GPa. c,d Predicted formation enthalpy of MgZn304, MgZnOs, and Mg3Zn04 with respect to 
decomposition into Bl-MgO and (c)B4-ZnO or (d) B3-ZnO below 8 GPa. 



unstable phase mixing^'^^, as reported by Sanati et aU^. The same 
conclusions are made by Seko et al.^^ and Malashevich et aV^, who 
both get that the formation energy of MgxZni_xO is positive. Are 
there any stable structures of MgxZni_xO alloy when ZnO and MgO 
adopt the same structure (51) at pressures greater than 9 GPa? 
However, for the ordered structures of Mg^Zni-xO under high pres- 
sure, few theoretical studies have been reported to the best of our 
knowledge. Therefore, a first-principles calculation is very interest- 
ing to investigate the high pressure stable structures of MgxZni _ xO, 
because they are the basis of all properties in theoretical studies and 
can also provide valuable information for experimental synthesis. 

Here, we investigate the possible stability of MgO -ZnO alloys 
under high pressure using the first-principles calculations, and the 
structures are obtained from a recently developed evolutionary algo- 
rithm for the prediction of crystal structures^^'^°. We also analyse the 
band gaps of these alloys. 

Results 

In this work, we report a few stable ground state structures of 
MgxZni _ xO. The phase stabilities of MgO-ZnO systems are inves- 
tigated by calculating the formation enthalpy of various MgxZni _xO 
alloys at different pressures. The formation enthalpy of MgxZni _xO 
is calculated by using fractional representation MgxZni _xO (0 < x < 
1) with respect to the decomposition into MgO and ZnO, as 

AH(Mg^Zni_xO) =H(Mg^Zni_xO) 

- [xH5i(MgO) + (l -x)H53,s4,Bi(ZnO)], 

where the enthalpies H for MgxZni _xO are obtained for the most 
stable structures as searched at the desired pressures. For MgO, the 



known structure B\ at our studied pressure is considered. For ZnO, 
transition from 54 to 51 takes place at 8.8 GPa as suggested by our 
enthalpy calculation (see Supplementary Fig. SI), which is in good 
consistency with the experimental findings pt ~ 8.7 GPa^^, pt ~ 
9.1 GPa^^ ox pt ^ 10 GPa^^ and first-principles calculations pt = 
9.3 GPa^\ It can also be seen that the B3 phase is slightly higher in 
energy than the B4 phase of ZnO. So, the enthalpies H for ZnO below 
8.8 GPa, B4, B3 and Bl phases are all considered, for the reason of the 
lowest energy, lower energy and cubic structure, and Bl-MgO-like 
structure, respectively, and above 8.8 GPa, only Bl-ZnO is consid- 
ered for the lowest energy and B 1 -MgO -like structure. The formation 
enthalpies of MgO-ZnO alloys under conditions of high pressure are 
depicted in Fig. 1 and Fig. 2. From Fig. 1 we noticed that MgZn304, 
MgZn02, and Mg3Zn04 made from Bl-MgO and Bl-ZnO are ther- 
modynamically stable at zero pressure because of their negative 
formation enthalpy, which is consistent with the results reported 
by Liu et aU^ While in the case of B4-ZnO (or B3-ZnO) and Bl- 
MgO, the formation enthalpies of all MgxZni _xO alloys we got are 
positive below 6.5 GPa (or 5.9 GPa), indicating the tendency for 
segregation of the ZnO and MgO, in agreement with the calculated 
results at ambient pressure of previous theoretical studies^^'^^'^°. 

Further detailed calculations above 8.8 GPa found a number of 
MgxZni _xO stable against decomposition into the constituent oxides 
(Bl-ZnO, Bl-MgO, or MgxZni _xO). The formation enthalpies are 
shown in Fig. 2 for pressures of 10, 40, 60, and 80 GPa, and detailed 
enthalpy calculations for the most stable structures allowed us to 
construct the P-x phase diagram of the MgO-ZnO alloys (see 
Fig. 3). When the pressure reaches 10 GPa or higher, the formation 
enthalpies of all MgxZni _xO alloys we predicted become negative. 
This is interesting trend for the number of MgxZni _xO alloy to 
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Figure 2 | Thermal stability of MgO-ZnO alloys. Predicted formation enthalpy of MgO-ZnO alloys with respect to decomposition into Bl-MgO and Bl- 
ZnO at (a) 10 GPa, (b) 40 GPa, (c) 60 GPa, and (d) 80 GPa, respectively. 



increase with increasing pressure. In this direction, high pressure is 
obviously the right tool to widen the range of compositions in which 
the stable structures can be obtained. To the contrary, the structures 
of Mg2Zn03 and Mg3Zn04 become metastable when pressure is 
increased above about 50 GPa. This resulted in a confirmation of 
the following structures as DFT predicted ground states under dif- 
ferent pressures: MgZn304 (Pm-3m, I4/mmm), MgZn02 (P4/mmm, 
I4i/amd), MgaZnsOs (C2/m), MgsZnOj (Cmcm), Mg3Zn04 
(Pm-3m, lAlmmm), Mg4Zn05 (J4/m), MgyZnsOg {Cllm), and 
MgeZnOy (i^-3). For all the predicted structures as shown in Fig. 4 



we computed phonon dispersions (see Supplementary Fig. S2) in the 
pressure range of 0-80 GPa, and found them to be dynamically 
stable. 

As can be seen from Fig. 2, at 10 GPa, the energetically "deepest" 
structures occur at x = 0.25, x = 0.5 and x = 0.75, which have been 
extensively studied by many authors^^"^^. Two structures were found 
respectively at the three compositions under pressure bellow 80 GPa, 
and the enthalpies difference are depicted in Fig. 5. It should be 
pointed out that the enthalpy differences are only within a few 
meV/formula. At 6.5 (or 5.9 GPa for B3-ZnO considered) to 
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Figure 3 | Pressure-composition phase diagram of the MgxZni _xO alloys. Red line represents the pressures at which MgZn304, MgZn02, and Mg3Zn04 
become stable against decomposition into B3-ZnO and other oxides (Mg^Zni-xO, or Bl-MgO). 
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Figure 4 | Crystal structures of the predicted MgxZni_xO alloys, (a) Pm- 

3m_MgZn304, (b) WmmmJsAgZn^O^y (c) I^il amdJsAgZnO 2y (d) P4/ 
mmm_MgZn02, (e) C/m_Mg3Zn205, (f) Cmcm_Mg2Zn03, (g) Pm- 
3m_Mg3Zn04, (h) I4/mmm_Mg3Zn04, (i) C2/m_Mg7Zn209, (j) 74/ 
m_Mg4Zn05, (k) P-3_Mg6Zn07. The large orange, middle gray and small 
red spheres represent Mg, Zn and O atoms, respectively. For more 
structural details, see the Supplementary Table SI online. 



21 GPa, Mg3Zn04 is stable in the Pm-3m structure (Ll2-type, as 
reported in Ref. 25), and then transforms into a tetragonal 14/ 
mmm structure (D022-type), which is the same as that reported in 
Refs.^^'^^. For MgZn304, the cubic Pm-3m and tetragonal 14/ mmm 
phases are found to be the most stable in pressure ranges of 8.5 (or 
7.9 GPa for B3-ZnO considered) -68 and 68-80 GPa, respectively. 
The Pm-3m MgZn304 (Ll2-type) is the same as that reported in 
p^g£g24,25 ]y[g2n02 is stable in the P4/mmm structure (same as that 
in Ref. 25) at 7.7 GPa (or 7.1 GPa for B3-ZnO considered), and 
transforms into a tetragonal MJamd structure (same as "40"-type 
reported in Ref.24) at 38 GPa, which is stable up to at least 80 GPa. 
Fig. 5(b) shows that both of the MgZn02 structures are more stable 
rather than in chalcopyrite structure^^'^^. It is emphasized that the 
structural researches of Refs^^'^^'^^ are all performed at atmospheric 
pressure. 

Detailed calculations and analysis of the other structures under 
higher pressure are performed and show other five stable MgO-rich 
MgO-ZnO alloys made from isostructural components (51-MgO 
and 51-ZnO): Mg2Zn03 (Fig. 4(f); space group Cmcm), stable in a 
pressure range of 28-53 GPa; and Mg7Zn209 (Fig. 4(i); space group 
C2/m), Mg4Zn05 (Fig. 4(j); space group J4/m), Mg3Zn205 (Fig. 4(e); 
space group C2/m), and MgeZnOy (Fig. 4(k); space group R-3), 
which are stable above 23 GPa, 41 GPa, 48 GPa, and 60 GPa, 
respectively. Among the five phases of MgxZni_xO alloy, the C2/m 
Mg7Zn209 is the same as Rl-Mg7Zn209 (C2/m) reported in Ref.^^. 
For the MgO-rich side of the phase diagram, we predict a stable 
hexagonal R-3 Mg6Zn07, which confirms the possibility that 
MgxZni_xO crystallizes in a hexagonal structure at high concentra- 
tions X, as recently reported by Shimada et aU^. This is a very inter- 
esting case due to the point that MgxZni_xO alloy with high 
concentrations x have long been thought to adopt rocksalt (cubic) 
structures. Especially, the MgO-ZnO alloys more should be 51 -like 
structure since both MgO and ZnO are stable in 51 phase above 
8.8 GPa. In addition, our calculated formation enthalpy data of 
metastable C2/m-Mg4Zn307 is slightly above the convex hull up to 
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Figure 5 | Calculated enthalpies as the function of pressure, (a) Enthalpy curves (relative to lA/mmm phase) for Pm-3m MgZn304, (b) Enthalpy curves 
(relative to 14 J amd phase) for P4lmmm MgZn02, (c) Enthalpy curves (relative to Immm phase) for Cmcm Mg2Zn03, and (d) Enthalpy curves (relative to 
141 mmm phase) for Pm-3m Mg3Zn04 by CASTEP code. The values in brackets represent the pressures at which MgZn304, MgZn02, and Mg3Zn04 
become stable against decomposition into B3-ZnO and other oxides (Mg^Zui-xO, or Bl-MgO). 
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Figure 6 | Band-gap variation of MgxZni_xO alloys as a function of MgO concentration at different pressures, (a) Calculated LDA band gaps, (b) 
Comparison of LDA band gaps with HSE06 values. 



80 GPa, which is different from their report on R2-Mg4Zn307 (C2/ 
m) as one of ground- state structures'^. 

Figure 6 illustrates the calculated band gaps as a function of 
Mg composition using the the local density approximation 
(LDA-CAPZ)^''^^ functional and screened hybrid functional of 
Heyd, Scuseria, and Ernzerhof (HSE06)^^'^^ for the ground-state 
MgxZni_xO structures under high pressure. From Fig. 6(a), it can 
be clearly seen that the gap increases with increase in the Mg content 
at the same pressure and increases with pressure increased at the 
same concentration. In addition, two different slops occur at about 
X = 0.75 of Eg-x curves for different pressures, which agrees with 
the findings reported in previous works^^'^^"^^, in which the change 
of slope near one composition value has been attributed to struc- 
tural phase from wurtzite to cubic with Mg content increased. 
Our calculated results show that no phase transition from 
wurzite to cubic or vice versa takes place at the point where the 
slope changes, and we consider that the band gaps of MgxZni_xO 
strongly dependent on the Mg content. The HSE06 functional results 
are depicted in Fig. 6(b), which show the same trend as the LDA 
values, with the difference of 0.5- LO eV higher value at the same 
pressure. 

Discussion 

We have performed a systematic theoretical investigation of the 
compositions and structures of the MgO-ZnO alloys and identified 
the stable ground- state structures as well as a number of metastable 
structures at 0-80 GPa over a wide range of Mg concentrations. Our 
results show that stable ordered ground- state MgxZni_xO alloys 
only occur at high pressure. As well as recognizing the previous 
theoretically observations for the MgZn304 (Pm-3m), MgZn02 
(P4/mmm, I4i/amd), Mg3Zn04 (Pm-3m, I4/mmm), and 
Mg7Zn209 (C2/m) at atmospheric pressure as high pressure stable 
phases, the new stable alloys MgZn304 {141 mmm), Mg3Zn205 (C2/ 
m), Mg2Zn03 (Cmcm), Mg4Zn05 (J4/m), and MggZnOy (R-3), have 
been predicted in the ground state under different pressure ranges. 
Among these structures, only MgeZnOy with a high Mg content 
crystallizes in the hexagonal structure under high pressure, and for 
the entire alloy composition of MgxZni_xO, there is no phase trans- 
ition from wurzite to cubic or vice versa takes place. Moreover, the 
band gaps of MgO-ZnO alloys with high MgO content increase more 
rapidly than low MgO content. 



The high pressure ground- state and metastable structures are all 
metastable at ambient pressure and may be synthesized under certain 
conditions (e.g., at high temperature). 

Methods 

Crystal structure prediction. The ground state structures of MgxZni_xO below 
80 GPa are predicted by the ab initio evolutionary algorithm USPEX'^'^'^". We have 
performed variable-composition simulation through the USPEX code in the range of 
0-80 GPa for (MgO)m-(ZnO)n system (n + m < 30), and determined the 
compositions of the MgxZni_xO. To confirm this and to obtain the most accurate 
results, we then renewed our structure search of stoichiometric MgmZnnOm + n with 
m and n values fixed by the USPEX code. 

Total energy calculations. The structural optimizations, electronic structure, 
phonon dispersion, and energy calculations for selected structures are performed 
with CASTEP code^^. The local density approximation (LDA-CAPZ)^^'^^ approaches 
of exchange-correlation functional is employed. The norm- conserving 
^seudopotentials*^" with cutoff energy of 700 eV are used, and a k-mesh of 0.03 X 2n 
A"^ within the Monkhorst-Pack scheme for sample the Brillouin zone, which ensures 
the error bars of total energies are less than 1 meV/atom. The hybrid functional 
HSE06^^'^^ implemented in the CASTEP code, together with norm-conserving 
pseudopotentials, and a cutoff energy of 600 eV, are used to calculate the band gaps of 
the stable structures at 10 GPa and 80 GPa. 
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